Winding Numbers, Complex Currents, and Non-Hermitian Localization 



OO 
c3 



Nadav M. Shnerb and David R. Nelson 
Lyman Laboratory of Physics, Harvard University, Cambridge, MA 02138 

(February 1, 2008) 

The nature of extended states in disordered tight binding models with a constant imaginary 
vector potential is explored. Such models, relevant to vortex physics in superconductors and to 
population biology, exhibit a derealization transition and a band of extended states even for a one 
dimensional ring. Using an analysis of eigenvalue trajectories in the complex plane, we demon- 
strate that each delocalized state is characterized by an (integer) winding number, and evaluate the 
associated complex current. Winding numbers in higher dimensions are also discussed. 
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There is a growing interest in the spectra of random 
non-Hermitian matrices [EJ. Considerable attention has 
focused on a particularly simple class of tight binding An- 
derson models with a constant imaginary vector poten- 
tial, inspired by the physics of vortex matter [0 . These 
models exhibit a sharp derealization transition even in 
one and two dimensions. Similar operators, represented 
by large real, asymmetric sparse matrices, arise in the- 
ories of population biology in random media with con- 
vection ||, and in many other contexts M. Delocalized 
eigenmodes arise in response to a sufficiently large asym- 
metry parameter, accompanied by eigenvalues which es- 
cape in conjugate pairs from the real axis into the com- 
plex plane [gj. 

Although the initial work relied heavily on numerical 
investigations || , there has been considerable recent ana- 
lytic progress, particularly in one dimension. Brouwer et. 
al. calculate explicitly finite size effects and the bubble of 
complex eigenvalues representing extended states in the 
center of the band in the limit of weak disorder || . Brezin 
and Zee Q have obtained exact results for the density 
of states for the special case of Lorentzian site random- 
ness embodied in the Lloyd model [Q. Goldsheid and 
Khoruzhenko discuss this model and show more generally 
that the eigenvalues are distributed along curves in the 
complex plane, providing analytic formulas relating the 
spectrum to the properties of a reference Hamiltonian 
with no asymmetry |t]]. The nontrivial behavior which 
results for large asymmetry parameters in two and three 
dimensions has been studied analytically via renormal- 
ization group calculations and a mapping onto Burgers 
equation ||. 

In this paper we study the eigenfunctions and complex 
currents associated with the band of extended states in 
one dimension. Unlike delocalized states in Hcrmitian 
disordered systems (where the eigenfunctions can always 
be chosen to be real) , we show that these complex eigen- 
functions are characterized by a conserved winding num- 
ber n even when the disorder is strong. Such topological 
quantum numbers can be used to label the eigenvalue 
spectrum e n (g), where g is the asymmetry parameter, 
i.e., the imaginary vector potential. A study of the eigen- 
value trajectories as a function of g then leads to complex 



currents J n = —i^§^, which determine the average tilt 
dr/dz of the vortex trajectory r(z) in a superconductor 
with columnar pins in the presence of a perpendicular 
magnetic field proportional to g [0. Assuming for sim- 
plicity a periodic box of size L z in the column direction, 
one has, for example, 



< 



dr(z) 
dz 



>= 



iLz/T 



-e n L z /T 



(1) 



where the brackets denote a thermal average at temper- 
ature T. We shall also argue that winding numbers can 
be used to label extended states for large asymmetry pa- 
rameters in two or more dimensions. 

The one dimensional non-Hermitian tight binding 
model "Hamiltonian" in a basis of N sites localized at 
positions {rj} reads, 



Hij = -Wj (e^'tfj+ij- + e 



(2) 



where we assume periodic boundary conditions and the 
{hj} are real asymmetry parameters. When exponen- 
tiated, this operator describes the transfer matrix for a 
flux line with columnar defects in a cylindrical shell [||. 
The Liouville operator describing population growth on 
a lattice in an inhomogeneous environment with convec- 
tion is given by C = —Ti. ||. The zero mean random 
potential Vj is chosen independently for each site, and 
arises from the variations of columnar pin diameters (vor- 
tex matter) or inhomogeneous growth rates (population 
biology). As discussed in the randomness in the off- 
diagonal hopping is due to the irregular spacing between 
columnar defects in superconductors. If Oj is the spacing 
between nearest neighbor columns j and j + 1, then the 
Wj = wq cxp(— 7£tj) and hj = gaj/a, where 7 and wq > 
are constants, a is the average lattice constant, and g 
is proportional to the field component which tries to tip 
vortices away from the columns. In models of population 
growth, Wj describes diffusion between inhomogeneously 
distributed population centers and hj represents a fluctu- 
ating convection velocity with average value g and second 
moment a. We impose periodic boundary conditions in 
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the space-like direction i.e., the N + 1-st site is identified 
with the first site of the sample. 

Upon carrying out a similarity transformation on the 
Liouville operator, C — SLS" 1 , where 



3-1 / • _ n 
Si,3 = &i,3 ex P[X! hk AT ^2 hk 



fc=l 



(3) 



fc=i 



one finds that the spectrum will be the same as that 
of a Liouvillian with a uniform asymmetry parameter, 
namely £' = —TL' with 
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FIG. 1. Eigenvalue trajectories in the complex plane. In- 
set: Angles entering Eq. (6) for a given e(g). 
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where g = Y^,j=i hj- Eqs. (|| - 1|) generalizes the gauge 
transformation used to describe localized states in, e.g., 
Ref. Q and || . The sample to sample fluctuations of g 
about its mean fall off like a/VN ||. As N grows one 
may thus replace the fluctuating quantity hi by a disor- 
der independent average value g, as we shall do in the 
rest of this paper. This result implies universality in the 
response to random convection in ID growth models - 
nonuniform convection velocities may be mapped into a 
uniform average velocity via the transformation (|s|) , pro- 
vided one applies this similarity transformation to the 
eigenfunctions as well. 

For a ID ring with random parameters {t^}, {u)i\ and 
g = 0, all the eigenfunctions of (Q) are real and localized, 
and its eigenvalues are real and discrete || . We assume 
a large but finite chain, such that although the spectrum 
is discrete, the length of the chain is much larger than 
the maximal localization length of an eigenmode. 

In a typical symmetric one-dimensional disordered sys- 
tem with g = 0, the localization length £ is largesy at 
the center of the band and smaller at the tails. The 
criterion for delocalization of the asymmetric system is 
Ka < g, where k = l/£ fl3|| . As a result, pairs of 
complex energies representing delocalized states first ap- 
pear as a "bubble" at the center of the band, which 
then spreads into the band tails. To study the com- 
plex currents associated with delocalized states we fol- 
low Refs. [5-7] and exploit the relation between the com- 
plex spectrum of the asymmetric problem with j / 
and the real eigenvalues of a "background" localized 
problem with g — 0, and same realization of disorder. 



The condition for a complex number e = e^ + ie/ to be 
an eigenvalue of the matrix C'(g) is Det[e — C(g)_ = 0. 
Eq. (Q) implies that e is an eigenvalue only if ]^|l2| 



Detle-C'(g = 0)] = l[(e-et) 



i=i 



N 



2[cosh( 3 A0 - 1] Yl 

3 = 1 



(5) 



where the {e"} are the (real) eigenvalues of the back- 
ground matrix C'(g = 0). To extract winding numbers, 
we first observe that the right hand side of (||) is real 
and positive for even N, while for odd N it is real and 
negative. As a result, each complex e should satisfy (see 
inset to Fig. 1) 



p-K, 



(6) 



with p = 2n for even N, p — 2n + 1 for odd N, where n 
is an integer, and the function cot^ 1 (x) varies from ir to 
as 1 goes between — oo and oo. 

As e; ^ 0, each term in Eq. (||) gives ir for every 
eigenvalue e° to the right of cr, and zero for each eigen- 
value to the left. To satisfy @ for even (odd) N and a 
given n, the eigenvalue must leave the real axis and en- 
ter the complex plane at the gap between the 2n-th and 
the 2n+l-th (2n+l-th and the 2n-th) eigenvalues of the 
g = "background" system. We call n the index of the 
trajectory en(g) + iei(g) in the complex plane. 

For odd N, we see immediately from (0) that the 
rightmost eigenvalue (with p — n — 0) must remain 
real, consistent with Perron- Frobenius theorem [n]]. The 
corresponding nodeless eigenfunction corresponds to the 
ground state of TL' = —C. For N even, particle- hole 
symmetry implies that both the rightmost and leftmost 
eigenvalues are always real. More generally, for a fixed 
value of n, the set of all [en{g) 1 £i(g)] satisfying Eq. (||) 
defines a curve in the complex plane, as illustrated in Fig. 
1. Henceforth, we assume N even for simplicity. 
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A more complete description of the eigenvalue trajec- 
tory results from taking the logarithm of the modulus 
of Eq. (|J). In the limit Ng >> 1 one finds a second 
constraint on e(g), namely, [5-7] 

M = -ln(f) + ^£ln(|6- e °|) (7) 

i 

where |e-e?| = ^/(e R - eV) 2 + ej and w = 2[E[f =1 2§L]*. 

This constraint is described graphically in Fig. 1, 
which shows schematically the level curves defined by 
(7) near the band center for three values of g. These 
are lines of constant potential for an equivalent 2d elec- 
trostatic problem with charges at the positions of the 
localized point spectrum for g = 0. When g is small, the 
constraint is solved by eigenvalue pairs on the real axis 
in the gaps between neighboring e°'s indexed by p = 2n. 
The departure of the eigenvalues from their g — values 
is initially exponentially small in the system size pj . As 
g increases, however, successive pairs of eigenvalues even- 
tually merge at a saddle point in the potential contours 
and detach from the real axis at right angle. For a given 
real energy ej?,, Thouless has defined an energy depen- 
dent inverse localization length k,(£r) for the associated 
Hermitian problem Jl2| , 

K ( e *) = -ln(|) + l]Tln(| ei? ,-e°|). (8) 

i 

Upon comparing with Eq. (m), we see that e/ becomes 
nonzero whenever \g\ > K(e R ) where (e'^,0) is the de- 
tachment point of the eigenvalue pair. 

As g increases above n{e' R ), Eqs. (JsJ) and ^ define a 
unique pair of complex eigenvalue trajectories e n (g) and 
e*(<?) for every value of n. Upon passing to the limit 
N — > oo , the spectrum {e°} for g — closes up, and 
is described by a density of states po(X). Eq. (Q) and 
(0) may then recombined into a single complex equation, 
namely 

J™ dX p (X) ln[e n (g) - A] = In (| e l*l) + (^j . 

(9) 

where la[e(g)-X] = ln[\e R +iei-X\]+icot~ 1 [(e R -X)/e I ]. 
In the limit N >> 1 and for a density of states func- 
tion symmetric around the special detachment point with 
e' R = 0, there is a purely imaginary trajectory of the form 
e(g) = iei(g) with n = N/2. For fixed g, Eq. (|) thus 
leads to an implicit formula for the "height" e™ ax of the 
bubble of complex eigenvalues in the center of the band, 
namely, § / dX p{X) ln[A 2 + ( e ™^) 2 ] = \g\ + m (f). This 
integral vanished, as expected [6], in the "one way" limit, 
g — > oo with ^e' 9 ' = 1. For other detachment points, 
the eigenvalue trajectories curve to the left or right as 
required by the constraint (6) (See Fig. 1). 

The analysis above suggests that the imaginary parts 
of most {e n (g)} diverge as \g\ — > oo. It is then expected 



that all eigenfunctions 4> n {j) are approximately plane 
waves, 4> n (j) ~ exp(ik n rj/a), with free particle eigen- 
value energies e n (g) = wcos(k n + ig) [|]||. For large 
\ei{g)\, Eq. (|) leads to 

w 

e n (g)^ -eM\g\ + ^n/N). (10) 

Comparison with the free particle spectrum at large \g\ 
shows immediately that the index n of the eigenvalue 
trajectory and the wavevector are related, k n = 2nn/N. 
These wave eigenfunctions spiral around the origin in the 
complex plane as one moves along the Id lattice of the 
tight binding model sites, leading to a well defined wind- 
ing number n. The winding number associated with the 
eigenvalue e* (g) in the lower half plane is then —n. 

As g decreases, the associated delocalized wave func- 
tion 4> n (j]g) must remain nonzero at every site: If 
<p n {j]g) were in fact exactly zero at some site i, it can be 
shown that the state is localized with a real eigenvalue 
and eigenfunction by mapping all effects of the asym- 
metry onto the special site i via a transformation like 
(^). Thus, the winding number must be preserved as g 
decreases, i.e., the winding number is a topological in- 
variant along an eigenvalue trajectory. The projection of 
such a delocalized eigenfunction is illustrated in the inset 
to Fig. 2. 

As N — + oo, we can replace the winding index n 
by the continuous variable k n — ^-n. Eq. (||) then 
shows quite generally that the complex spectrum is a 
function only of the combination g + ik, i.e., e n (g) = 
t{g + ik). It then follows from the Cauchy-Riemann 
relations that all eigenvalue trajectories are at right 
angles to the lines of constant g, de j 9 ^ = —i de }^ ■ 
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FIG. 2. Spectrum and winding numbers for N — 1000 
with Vj uniformly distributed in the interval [—1, 1], and 
g — 0.4. Inset: Projection onto the complex plane of 
eigenfunction with N = 200, g = 1, and n = 100. 

An explicit formula for the complex current results, 
moreover, from differentiating Eq. (||) with respect to g, 
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J[ e n(g)} = -i 



dg 



d\ 



Po(A) 



e n (g) - A 



= -iG^ l [e n {g)] (11) 



Evidently, the current associated with a particular 
complex eigenvalue e n {g) is determined by the Green's 
function of the g = problem. As an application of 
our results, we note that the disorder averaged Green's 
function for the Lloyd model, with Lorentzian site dis- 
order, Pr^) = Ijy^ is @@ < G (E) >= 

— 1/2 

[{E + 17) — t] , and complex energy spectrum for ex- 
tended states with Ime„ > 0, e n (g) = w cos(kn + ig) — 17. 
It then follows immediately from Eq. (pi]) that the 
disordered averaged inverse complex current J n = | < 
J- 1 >- x I is 



Jn = —i w sm{k n + ig) 



(12) 



in agreement with a direct calculation of the current from 
the spectrum itself. If we focus attention on Im J as a 
function of Re e, we obtain ImJ = — w tanh(g) Ree, in 
excellent qualitative agreement with the results obtained 
for the imaginary current with bounded disorder in Ref. 
Pj for different values of g. Although we have checked 
our results with the Lloyd model, we stress that Eq. (|Tl| ) 
expressing the complex current in terms of the Green's 
function of the Hermitian reference system is much more 
general. 

Fig. (2) shows a typical eigenvalue spectrum for £ in 
one dimension superimposed on the winding numbers as- 
sociated with the extended eigenfunctions. The winding 
numbers are exactly ±N/2 at the band center, and their 
magnitudes decrease monotonically as one moves toward 
the upper edge of the band, i.e., toward the lowest en- 
ergies of the corresponding Hamiltonian. The winding 
numbers remain finite up to the mobility edge separat- 
ing complex and real eigenvalues, but become undefined 
in the band tails, where all wavefunctions are real. In 
this case, the winding index classification can be replaced 
by counting the nodes in the localized wavefunction [fl2[ . 
The size of the inner excluded region of the projected 
wavefunction shown in the inset shrinks to zero at the 
mobility edge. The winding number at the mobility edge 
decreases smoothly to zero as g increases and all states 
become delocalized. 

Given that a well defined topological quantum num- 
ber (which plays the role like that of the momentum) 
exists for the extended states in the one dimensional non- 
Hermitian localization problem, it is interesting to specu- 
late whether similar topological invariants exist for delo- 
calized states in higher dimensions. For large g in D > 2, 
disorder produces nontrivial changes in the spectrum, but 
the eigenfunctions are still slightly perturbed plane waves 
H . For NxN lattices with toroidal boundary conditions, 
it seems likely that the winding number (n x , n y ) embod- 
ied in the wavevector k = (k x ,k y ) — (2im x /N,2im y /N) 



is initially preserved as g decreases from large values. 
However, no detailed analysis exists of how these ex- 
tended states in D = 2 become localized as their eigen- 
values approach the real axis, and we cannot be certain 
if the winding number remains strictly invariant along a 
trajectory. In this context, it is interesting to note that a 
recent analysis of delocalized wavefunctions for the non- 
Hermitian problem in D = 2 predicts algebraic decay of 
the wavefunction amplitude |l4| , similar to the quasi-long 
range order below the vortex unbinding transition in sta- 
tistical mechanics. Vortex pairs may also be involved in 
changes of the vector winding number of delocalized wave 
functions for d = 2, similar to the decay of supercurrents 
in Helium films at finite temperature Jig] , 
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